Electron-Doped Manganese Perovskites: The Polaronic State 
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Using the Lanczos method in linear chains we study the ground state of the double exchange 
model including an antiferromagnetic super-exchange in the low concentration limit. We find that 
this ground state is always inhomogeneous, containig ferromagnetic polarons. The extention of the 
polaron spin distortion, the dispersion relation and their trapping by impurities, are studied for 
diferent values of the super exchange interaction and magnetic field. We also find repulsive polaron 
polaron interaction. 
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The discovery of 'colossal ' magnetoresistance (CMR) 
(jjj together with its many unusual properties has received 
considerable attention lately. 

Experiments have revealed very rich phase diagram in- 
terpreted in terms of magnetic ferro, antiferro, canted, 
polaronic, and non-saturated phases . Charge ordered 
phases also have been found ||. The phase diagram, 
as a function of concentration x, temperature, magnetic 
field, or magnitude of the superexchange interaction is 
not quite clear yet for the different compounds. The 
metallic phase can be reached by hole doping of the par- 
ent compound LaMn03, by substituting La for divalent 
alkalies, Pb, or by stoichiometry changes. Very recently, 
neutron scattering experiments have been interpreted in 
terms of polaronic droplets || . Much less is known about 
the electron doped compounds where doping does not 
seem to produce metallization. 

From the theoretical point of view, the pioneering work 
of de Gennes H proposed a canted phase to resolve the 
competition between the ferromagnetic double interac- 
tion introduced by the presence of itinerant holes and 
the superexchange interaction. Recently, several contri- 
butions to this problem have been reported. Arovas and 
Guinea p], studied this problem using a Schwinger bo- 
son formalism to obtain a phase diagram showing several 
homogeneous phases and pointing out that phase separa- 
tion replaces the canted phase in a large region. Indeed 
phase separation appears in several numerical treatments 
of the problem ||. In other analytical treatments, more 
adequate to treat local instabilities, non-saturated local 
magnetization states have appeared at zero temperature 
0. M.Yu.Kagan et. al. || have studied the stability of 
the canted phases against the formation of large ferro- 
magnetic 'droplets' containing several particles and they 
conclude that the formation of droplets is favored in the 
ground state. The variety of results obtained from the 
different approaches points to the need of clarifying the 
picture and testing the results. 

In this work, for the first time, we find the low en- 
ergy quasiparticles and characterize their structure and 
dispersion relation in the low concentration limit. These 
quasiparticles correspond to the electron followed by a 
ferromagnetic local distortion (ferromagnetic polaron) in 
the antiferromagnetic (AF) background. The dispersion 



relation is dominated by k — > k + ir scattering due to the 
presence of AF order. In order to make a conection with 
transport properties, we also study the the tendency to 
localization of these polarons in the presence of impuri- 
ties and magnetic field. 

To render evident the nature of the ground state, we 
resort to the Lanczos method, which is free from approxi- 
mations. The Hamiltonian is simplified to a single orbital 
per site, no lattice effects are considered, and we have to 
reduce to one dimensional chains. However our results 
provide a simple picture that, we presume, can be put 
to test of the dilute limit of electron doped systems. In 
these systems, the limitations of the model Hamiltonian 
may not be as stringent as in the hole doped systems 
for the following reasons: the lattice structure is more 
symmetric so Jahn Teller distortions should play a less 
important role, the large in-site Coulomb repulsion in- 
hibits double occupation so that it may be possible to 
describe the physics by the use of a single effective or- 
bital, and finally the antiferromagnetic structure of two 
interpenetrating lattices can be properly described in one 
dimension. 

In order to describe the manganites we consider two 
degrees of freedom: localized spins that represent the ti g 
electrons at the Mn sites, and itinerant electrons that hop 
from e2 g Mn orbitals to nearest neighbor &i g orbitals. 
The model Hamiltonian includes exchange (J) energies, 
an antiferromagnetic interaction between localized spins 
(K) and a hopping term of strength t which we will use 
as energy unit hereafter. It reads: 



H = -J h Y,S i -a i +K J2 S '' S J 

i <i,j> 

*« (c£ • c. ja + h.c.) 



(1) 
(2) 



<i,j>,cr 



where n* 



c~l a Ci<j , and cf a , Ci a creates and destroys 



an itinerant electron with spin a at site i, respectively. 
Si and <Ji are the localized and itinerant spin operators at 
site i, respectively. In what follows we take large values 
of J/j which prevents double occupation and makes the 
on-site Coulomb interaction U irrelevant. This model 
has been studied numerically for finite concentration in 
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reference 0] , and in the absence of AF coupling (K = 0) 
in reference Q . In this paper we focus on the dilute limit. 




FIG. 1. We show the correlation function < noSjSj+i > 
for a 16 sites chain, with A = 10 and different values of K. 
The maximun value of S for each value of K is indicated in 
the figure. One can observe that the total S and the exten- 
sion of the magnetic distorsion increase as K decreases. The 
oscillatory behavior found at K = 1 is a consequence of the 
weakening of the antiferromagnetic links around the charge. 
The inset shows < SjSj+i > for a 16 sites Heisenberg chain 
where the link at site zero is a factor 2 smaller than the rest. 

We first investigate the homogeneity of the solutions 
for different values of the antiferromagnetic interaction K 
(in units of t). To this end, we calculate the ground state 
with one electron added for chains of different sizes up 
to N = 20. With the aim of looking for spin distortions 
around the charge, we calculate a correlation function 
which makes such a situation evident : < riiSjSj+i >. 
Because of translational symmetry the results depend 
only on \i — j\. The results are shown in Fig. 1 where 
we plot N < noSjSj+i > vs j, where TV is the number of 
sites. As it can be seen in Fig. 1, for large j, this correla- 
tion function takes a value very close to the one obtained 
from the Bethe ansatz solution for the Heisenberg chain, 
< SjSj+i >= 0.443. The extension and the magnitude 
of the spin distortion around the particle increases as K 
decreases. The oscillations observed in the curve corre- 
sponding to If = 1 are also observed for larger values 
of K . They are a consequence of the weakening of the 
antiferromagnetic links around the charge position which 
produce a sort of local spin dimerization. To prove this 
point, we show in the inset the nearest neighbors spin- 
spin correlation functions for a Heisenberg chain of the 
same size where the link at site zero is a factor of two 
smaller than the rest. 

However as K decreases, it is difficult to find an ade- 
quate approximation to describe the large polaronic dis- 
tortion. In order to obtain the effective mass of these po- 
larons, we investigate the dispersion relation for charge 



excitations. To this end we calculate the lowest energy 
state for different values of the momentum k = 2irn/N 
within the subspace where the total spin is that of the 
ground state. In Fig. 2 we show the dispersion relation 
scaled to the thermodynamic limit for K = 3, Jh = 100; 
K = 1, J h = 10; and K = 0.3, J h = 10. 
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FIG. 2. Dispersion relation scaled to the thermodynamic 
limit for different sets of parameters. In full line we show 
a fit with the expression discussed in the text for the limit 
Jh » k > t with A = 1.4 ~ K/2. In dotted line fit with 
A = 1 and t = 0.75 for K = 1, and A = 0.6 and t = 0.23 for 
K = 0.3. 

We start analyzing the dynamics in the regime where 
(Jh » K >> t). In that case, the charge moves as a 
spin one (S). The effective hopping resulting from the 
projection of the hopping term onto the reduced S = 1 
Hilbert space is: tPij(£iSj + 1/2), where Py is the per- 
mutation operator between sites i and j. We can picture 
the movement of the particle, in this limit, as going from 
a state |t4fl"4T|t° an intermediate state | t I T 
t I, and finally to J. | 4 T 4 1T 4, where ft (0) repre- 
sents the S z = +1(0) components of the spin S = 1. 
Thus, in order to move, the charge has to hop to the 
nearest neighbor, via a spin flip process, through states 
that differ in energy by A « K/2. It can be easily veri- 
fied that the efective hopping of this process is equal to 
t e f —tj The dispersion relation given by this dynam- 
ics is: A/2 ± ^J(A/2) 2 + 4t 2 ef cos 2 (fc) . The expression 
corresponding to the lower band is plotted with full line 
in Fig. 2 and compared with the numerical result for 
K = 3 and Jh = 100. This expression is valid in general 
for a particle moving in an antiferromagnetic background 
where scattering between k and k + n states dominates 
the dynamics of the particle (dotted lines in Fig. 2). 

In the case where K >> Jh >> t the spin distortion 
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can be neglected and the particle propagates in an an- 
tiferromagnetic lattice. The Hund interaction alternates 
the site energy of the propagating particle so that the dif- 
ference between the two sublattices is given by A = Jh{< 
(jjSj+i > - < VjSj >) = J h (< SjSj+i > - < (jjSj >) 
where we approximate < <JjSj >« 1/4 its value at the 
triplet state, and < ajSj+i >~ < SjSj + i >= In 2 — 1/4 
, the Bethe ansatz value. Using these values we find 
A = 0.19J/J. In this case t e f is equal to t. 



show in full line the exponential fit of the different curves 
showing the change in localization length. 
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FIG. 3. Effective bandwidth as a function of K. The 
change of regime as K increases trough t is evident. 

When J h » t > K, the magnetic distortion around 
the charge is large and the effective hopping is dominated 
by the overlap between the magnetic distortions about 
the nearest neighbors sites. This last effect dominates the 
polaron effective mass. Therefore, the mass of polarons 
increases when K decreases, as obtained in Fig. 2, where 
t e f decreases from 0.75 for K = 1 to 0.23 for K — 0.3, in 
agreement with the above results showing that the spin 
distortion around the charge increases in magnitude and 
extension when K decreases. In Fig. 3, we calculate the 
bandwidth for several values of K. One can distinguish 
clearly two regimes: K « t and K » t, the first 
corresponds to a large magnetic distortion and the second 
corresponds to a smaller one according to Fig. 1. 

In order to test how robust is the polaronic description 
of the results, we pin the polaron to a site by lowering in 
eo the diagonal energy at site zero. This may be relevant 
to the real materials since the doping process necessarily 
introduces some disorder. Since we are treating a linear 
chain, this always localizes the particle, but the localiza- 
tion length should be very different for different effective 
masses, so that a small eo localizes much more the po- 
laron for low values of K than for larger ones. This is 
shown in Fig. 4 where we plot < rii > around site zero 
for different values of eo and K. For comparison we also 
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FIG. 4. Figure 4. Charge localization, we show < m > as a 
funcion of position in a 16 sites chain where the energy at site 
zero is lowered from the rest by 0.05 units of t. Full lines are 
exponential fits to the curves with the values of localization 
length A indicated in the figure. The short localization length 
of the lower K curves indicate effective bandwidths of the 
order of the energy change. 

In Fig 5, we show the change in the values of < rii > 
for different magnetic fields It can be seen that the lo- 
calization of the polaron decreases with magnetic field 
as a consequence of increased effective hopping between 
nearest neighbors. A fact that may be important for the 
transport properties of these systems since it implies a 
negative magnetoresistive behavior for conductivity due 
to hopping between localized states [H 
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FIG. 5. Figure 5. Effect of magnetic field on charge local- 
ization. We show < Hi > for different values of an external 
magnetic field for the same chain as in FIG. 4. 
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Finally, in Fig. 6 we present some of the results for the 
two particles ground states. The charge-charge correla- 
tion function, normalized to the non correlated (Jh = 0) 
case, clearly shows repulsion between the particles. This 
long range repulsion increases with the magnetic distor- 
tion indicating its magnetic origin. 
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FIG. 6. Figure 6. charge-charge correlation function for 
two particles in a ten site chain for A = 10, and different 
values of K. Values are normalized to the uncorretated values. 

The fact that polarons repel each other points to a pic- 
ture of the electron doped systems similar to that pro- 
posed originally by de Gennes Q of 'self-trapped elec- 
trons'. Further discussion on the non-diluted limit, will 
be postponed for a later publication where newer results 
will be shown. 

In summary, we have investigated the possibility of a 
non uniform ground state in a model Hamiltonian us- 
ing the Lanczos technique. The model describes chains 
of localized spins coupled antiferromagnetically on which 
electrons are added. These electrons suffer a strong fer- 
romagnetic interaction with the local spins and can hop 
from site to site. Assuming the model adequately de- 
scribes the physics of electron doped manganites, the 
results presented here point to a picture of these sys- 
tems where heavy one electron polarons dominate the 
magnetic and transport properties. Their masses depend 
strongly on the relation between the hopping energy and 
the AF superexchange interaction. Clearly, the doping it- 
self will localize the polarons so that transport will result 
from hopping between pinned sites. Negative magnetore- 
sistance should appear as a consequence of the decrease 
of the pinning energy with magnetic field JlO[ ] . 

Adding two particles to the chains we find long range 
repulsion between them. This long range repulsion could 
give raise to charge ordering. 

Finally, we would like to point out that the order of 
oxygen vacancies in CaMnO^-s makes real the possibil- 
ity of one dimensional electron paths in these materials 
|pr|. We hope that our results will stimulate more ex- 



perimental and theoretical investigations on the electron 
doped manganites. 
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